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ABSTRACT 



An algorithm is developed and evaluated for discriminating among clouds, snow 
cover and clear land. The multispectral technique uses daytime images of AVHRR 
channels 1 (0.63nm). 3 (3.7jim) and 4 (ll.Ojim). Reflectance is derived for channel 3 
by using the channel 4 emission temperature to estimate and remove the channel 3 
thermal emission. Separation of clouds from snow and land is based primarily on this 
derived channel 3 reflectance. Using this technique, observed reflectance in channel 3 
is 2 to 4 percent for snow, 3 to 10 percent for land, 2 to 27 percent for ice clouds and 8 
to 36 percent for liquid clouds. These values overlap for thin cirrus and snow, so the 
routine then attempts analysis of cirrus based on its different transmissive properties 
between channels 3 and 4. Six images were analyzed and the total cloud cover was 
verified against a total of 1 10 conventional surface observations using the standard 
categories of clear, scattered, broken and overcast. The routine was quite successful, 
with the analyzed sky cover being within category for 55 percent of the stations, one 
category different for 33 percent, 2 categories different for 9 percent and 3 categories 
different for 3 percent of the stations. A major remaining problem is discrimination 
between ice clouds and snow cover due to the great similarity of reflective properties of 
these two surfaces. 
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I. AUTOMATED CLOUD ANALYSIS 



A. INTRODUCTION 

One of the critical factors in the conduct of military operations is knowledge of 
the environment in which they are to take place, particularly regarding cloud cover. 
Clouds can dramatically affect such missions as aerial intelligence gathering, air 
refueling and tactical employment of weapons that use the visible-infrared portion of 
the electromagnetic spectrum to designate targets. Cloud cover analyses can be 
accomplished with synoptic weather observations but data-void areas severely limit 
their usefulness. Techniques have been developed to produce cloud analyses using 
satellite imagery, where human image analysts make subjective interpretations of cloud 
cover based on brightness and texture contrasts. Manual analysis takes a great deal of 
time to accomplish and quick response for large geographical areas in critical situations 
is very difficult. The Air Force Global Weather Central (AFGWC) developed an 
automated cloud analysis system in order to introduce objectivity to the analysis 
procedure based on sound physical principles, and also to address the Air Force's 
mission requirement for rapid cloud analyses. The 3-Dimensional Nephanalysis 
(3DXEPH) system became operational in January 1970. This pioneering effort did not 
have the benefit of experience from other systems. Various improvements have been 
made over the years, and in 1984 a second-generation model called the Real-Time 
Nephanalysis (RTNEPH) became operational. It is the only known automated system 
capable of processing and interpreting the tremendous volume of satellite data that are 
available today and integrating them with conventional information to provide a high- 
resolution. three-dimensional cloud analysis data base for the entire world. 

In addition to global analysis models, efforts have been made to develop 
automated analysis routines for minicomputers that can be used in a theater or 
battlefield environment. Wash et al. (1985) developed a cloud and precipitation 
analysis program for an interactive minicomputer system which uses geostationary 
infrared and visual data. This type of analysis package could be used in a weather 
station or a tactical van to provide cloud cover analyses and forecasts to operational 
commanders on a real-time basis. 
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B. SNQW/CLOUD DISCRIMINATION 

Various long-term problems have plagued automated systems, one of the most 
significant being discriminating low clouds from snow cover. Middle and high clouds 
can be detected against snow cover due to the thermal contrast detectable at infrared 
wavelengths, but in the regions of the electromagnetic spectrum used most often in 
satellite analysis, low clouds and snow cover have similar radiometric properties. In 
the visible portion of the spectrum both have high albedos, and in the infrared portion 
both have similar thermal properties. Thus there is little contrast in images at these 
spectral wavelengths. 

Bunting et al. (1977) discussed the reflectance properties of snow and clouds at 
visible and near-infrared wavelengths, using data from the Earth Resources Experiment 
Package which flew aboard Skylab. They analyzed the imagery in pairs: one set in the 
visible spectrum and the other in the near-infrared spectrum (1.55jtm to 1.75nm). The 
two sets showed high reflectance for snow, water clouds and ice clouds in the visible 
spectrum, but in the near-infrared the reflectance was high for water clouds, medium 
for ice clouds and very low for snow. They concluded that these differences could be 
the basis for an automated routine which makes an analysis decision based on the ratio 
of visible reflectance to near-infrared reflectance. Bunting and d'Entremont (1982) 
tested six automated classifiers on data from a special Defense Meteorological Satellite 
Program (DMSP) sensor in the wavelength band from 1.51^m to 1.63jim, along with 
visual and infrared data. They used the classifiers successfully to distinguish water 
clouds, ice clouds, snow cover and other cloud-free surfaces. Their technique was so 
successful that plans are now being made to add this sensor as an operational 
component of the DMSP satellite. 

Kidder and Wu (19S4) showed that there is a contrast in brightness temperature 
between snow and low clouds in daytime imagery of channel 3 (3.7nm) of the National 
Oceanic and Atmospheric Administration (NOAA) polar-orbiting Advanced Very High 
Resolution Radiometer (AVHRR). They related this contrast to solar reflection, which 
is a substantial component of the channel 3 radiance measurement along with thermal 
emission. Since low clouds and snow have similar emission temperatures, this contrast 
is the result of a difference in their solar reflection at 3.7nm which, in theory, is just a 
few percent for snow and about 20 percent for clouds. This results in a warmer 
brightness temperature for clouds (darker image) than for snow (whiter image). 
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Fye (1978) documented the 3DNEPH system and discussed the problem of 
snow cloud discrimination. Cloud cover is overestimated in areas of new snow and in 
areas where the ice analysis contains too little ice. Underestimation or no analysis of 
clouds results when snow melts rapidly or the ice analysis contains too much ice. 
Presently, the system uses the AFGWC Snow Cover Model and grid points that have 
snow or ice as a background are treated as missing for the cloud analysis routine. A 
solution to this problem has not yet been implemented on the RTXEPH system. This 
has a significant effect on the analysis in the late fall through early spring because 
about 30 to 35 percent of all the model grid points in the northern hemisphere are 
snow and or ice covered. 

C. THESIS OBJECTIVES 

The first objective of this thesis is to calculate and report solar reflectance in 
AVHRR channel 3 using a method that estimates the thermal emission from channel 4 
and removes this portion from the channel 3 radiance measurement. Derived values 
will be compared to theoretical 3.7jim reflectance for snow cover, land and clouds. 
During the 1990's both XOAA and the Department of Defense plan to launch 
satellites with a 1.6jim sensor for daytime snow/cloud discrimination, and the XOAA 
satellite will have a 3.7nm sensor for use at night. Even with this change, there are 
good reasons to pursue this work now. First, by the time the new sensor flies there 
will be over 10 years of archived daytime 3.7fim imagery. Second, since snow and 
cloud reflectances are similar at 1.6fim and 3.7jtm, a data set of refiectances from 
channel 3 would be very useful in testing software for the new sensor. 

The second objective of this thesis is to use calculated channel 3 reflectance to 
develop an automated cloud analysis routine that separates snow cover, clouds and 
snow-free land. It will be designed to handle various combinations of cloud and 
surface features, but the focus is on the separation of low clouds and snow. The 
routine will not rely on sophisticated statistical analysis or on artificial image 
enhancement techniques, but rather on basic physical principles of radiative transfer in 
the earth's atmosphere. In addition to similar visible reflective properties, low clouds 
and snow have nearly identical infrared thermal properties. Because of this and the 
results of the previously mentioned studies, the routine is based primarily on the 
reflective properties of these surfaces in AVHRR channel 3. A multispectral technique 
is developed and evaluated using visible and infrared data along with channel 3 to 
process daytime AVHRR images. 
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In Chapter II the necessary theoretical background is presented. The details of 
satellite data processing are presented along with a theoretical discussion of the 
reflection and emission characteristics of the cloud and surface features to be analyzed. 
A discussion of the satellite data and the techniques of data analysis are presented in 
Chapter III along with observed values of channel 1 and channel 3 reflectance for 
these surfaces. Finally, the analysis routine is developed using these measured 
reflectance values. Chapter IV discusses the performance of the analysis routine when 
applied to six images containing various combinations of snow and cloud cover, and 
sun-satellite geometry. Data were taken from two different NOAA satellites to 
illustrate the general utility of the routine and to provide additional statistical 
independence to the verification of the analyses. Total sky cover on the analyses is 
verified against the available surface observations using the standard conditions of 
clear, scattered, broken and overcast. The reported statistics include the percentage of 
stations where the analyzed sky cover matched the observed sky cover, and the 
percentage of stations with a one-category, a two-category and a three-category 
difference between the analyzed and observed sky cover. 
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II. THEORETICAL FOUNDATION 



A. DISCUSSION 

In general terms, detection of cloud cover on a satellite image depends on the 
contrast between the clouds and their background. At visible wavelengths this contrast 
is in terms of reflectance differences, and at infrared wavelengths this contrast is in 
terms of differences in brightness temperature. At middle infrared wavelengths 
(AVHRR channel 3) the contrast is a mixture of difference in reflectance and difference 
in brightness temperature. This chapter presents the theoretical foundation for the 
algorithm developed in this thesis. Included here is a theoretical discussion of the 
reflective and thermal characteristics of the features to be analyzed, the satellite data 
processing to derive these properties and a preliminary assessment concerning how the 
routine might work based on the theoretical characteristics of the features. See 
Appendix A for a list of symbols and constants used. 

Two visible and two thermal infrared channels are available on the AVHRR 
instrument. In addition to these, channel 3 is at a wavelength (3.7[im) that has 
contributions from both thermal emission and solar reflection. Channel 1 (0.63f.lm) is 
completely in the visible spectrum while channel 2 (0.87nm) extends into the near- 
infrared portion of the spectrum. Channel 1 is the visible channel used in this study. It 
will be shown later that the method of deriving the reflectance in channel 3 depends on 
channel 4 (1 l.O^im), so it is the infrared channel used in this study. 

Throughout this thesis, the satellite-measured radiance during daytime is 
approximated as follows: 

L = €B(T) + r(0 o ,e,<p)Icos0 o (2.1) 

The first term on the right hand side is the contribution to the measured radiance from 
thermal emission of the viewed surface assuming the transmissivity between it and the 
satellite is 1.0. The amount of radiance from thermal emission reaching the satellite is 
determined by the emissivity (£) of the viewed surface. The Planck function relates the 
emitted monochromatic intensity with the wavenumber and temperature of the 
emitting surface as follows: 



15 



B(v.T) - 



( 2 . 2 ) 



2hcV 

exp(hcv/KT) - 1 

The second term on the right hand side of Eq. 2.1 is the contribution to the measured 
radiance due to solar reflection. The amount of reflected solar radiance reaching the 
satellite is determined by the incident solar radiance (I) which is weighted by the cosine 
of the solar zenith angle (0 q ). and the reflectance (r) of the viewed surface. The 
directional reflectance is a function of the solar zenith angle, the satellite zenith angle 
(0) and the horizontal angle between them (q>). This sun-satellite geometry is 
illustrated in Fig. 2.1. Isotropic reflectance is related to the directional reflectance by 
the anisotropic reflectance factor (f). This relationship is expressed as (Taylor and 
Stowe, 19S4): 



r 3 (9 0 ,9.(p) 



(2.3) 



The anisotropic reflectance factor arises due to the directional dependence of the 
scattering mechanisms. It is the ratio between the radiant exitance in a given direction 
assuming the surface reflects isotropically, and the actual radiant exitance. For 
example, a radiance measurement in a given direction with a factor of 1.0 gives the 
correct radiant exitance, while a factor of 1.5 with the same measurement means the 
isotropic assumption would yield a 50 percent overestimate of radiant exitance. 

It is assumed that any viewed clouds are optically thick so that upwelling 
radiance from below them does not reach the satellite. This assumption breaks down 
when considering thin cirrus clouds and this will be discussed later. 

B. REFLECTANCE 

This section reviews the results of various theoretical studies concerning the 
reflective properties of the surfaces to be analyzed. These values then will be used for 
comparison to observed values in the data analysis section of Chapter III. 

1. Liquid Clouds 

The reflectance of liquid clouds is dependent upon cloud optical thickness and 
sun-satellite geometry. Tables 1 and 2 show the results of reflectance calculations from 
the model of Shettle and Weinman (1970), which is based on Mie theory and the delta- 
Eddington approximation. The values show the expected extremes of reflectance for 
clouds of various physical properties which are coded in the tables as follows: 
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